Density of states of colloidal glasses 
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Glasses are structurally liquid-like, but mechanically solid-like. Most attempts to understand glasses 
start from liquid state theory. Here we take the opposite point of view, and use concepts from solid 
state physics. We determine the vibrational modes of a colloidal glass experimentally, and find soft low- 
frequency modes that are very different in nature from the usual acoustic vibrations of ordinary solids. 
These modes extend over surprisingly large length scales. 



PACS numbers: 47.27.Cn, 47.20.Ft, 47.50.+d 

The glass transition is perhaps the greatest unsolved prob- 
lem in condensed matter physics [ 1 ] : the main question is how 
to reconcile the liquid-like structure with solid-like mechani- 
cal properties. In solids, structure and mechanics are related 
directly through the vibrational density of states (DOS) of the 
material [2]. It then seems important to obtain directly the 
density of states of a glass that is solid mechanically, but has 
no crystalline ordering. We do so for (colloidal) hard spheres 
that are known to undergo a glass transition upon increasing 
volume fraction. 

Recent theory shows that randomly packed hard spheres 
form a system with no redundant mechanical constraints. Be- 
cause of this, a small perturbation, like for example break- 
ing a particle contact, may induce rearrangements at all 
scales IH 0, IH : the system is 'critical' in this sense 10,0]. 
It has been argued that the marginality of the system has dra- 
matic consequences for the vibrational spectrum |8y: in stark 
contrast to well-known phonon modes in solids, which have 
been observed in previous studies fl [Tol[Tll[T2Tl . a broad band 
of floppy modes emerges in the vibrational density of states, 
with a gap at low frequencies that disappears as the pressure 
becomes infinite. 

In this Letter we measure, for the fist time, vibrational 
eigenstates of a colloidal glass/supercooled liquid system. We 
consider colloidal hard spheres that are subject to thermal ag- 
itation, allowing us to follow the 'vibrational' motion of par- 
ticles that are trapped in cages constituted of neighboring par- 
ticles. This is done at different random dense packing con- 
figurations around the glass transition. We obtain information 
on the DOS from a normal-mode analysis of particle displace- 
ments measured using confocal microscopy. We find that the 
vibrational spectrum has many soft low-frequency modes J^], 
more abundant and very different in nature from the usual 
acoustic vibrations of ordinary solids. This results in an 
anomalous low frequency peak in the density of states which 
approaches zero frequency as one goes deeper into the glass 
phase. The observed soft modes are collective 'swirling' par- 
ticle motions that extend over surprisingly long length scales. 

Hard-sphere colloidal systems exhibit a glass transition at 
which the dynamics becomes very slow around a volume frac- 
tion of <f> g iass = 0.58 [13]. We use poly-methylmethacrylate 
(PMMA) particles, of radius, a = 1.3/im, which are dyed 



with rhodamine and are sterically stabilised to prevent aggre- 
gation. The rhodamine dye makes it possible to follow the 
Brownian motion of all the individual particles on a plane us- 
ing a confocal microscope. There is a polydispersity of about 
5% in the particle size to prevent crystallization in the sys- 
tem. We use a mixture of cyclohexylbromide and decalin as 
our solvent, to match the density and index of refraction of the 
particles. The constituents and preparation procedure leads to 
particles that are indistinguishable from hard spheres. To ver- 
ify this, we have measured the crystallization density, which is 
a very sensitive measure for deviations from hard-sphere be- 
havior, and found it to agree to within a fraction of a percent 
to that of true hard spheres. 

Two-dimensional images were acquired in a field of view 
of 100 fim x 100 /im. The time interval between each im- 
age frame is 0.05 sec, which is approximately 1/10*' 1 of the 
Brownian timescale tb, tb = rjd 3 /ksT w 0.75sec, where T) 
is the solvent viscosity of the suspension. Following around 
2000 particles of radius 1.3/im in real time using a fast con- 
focal microscope (Zeiss LSM live) allows us to reconstruct all 
particle trajectories (Fig. Q]). 

The mean square displacement per particle < (Sr) 2 > as a 




FIG. 1: A two dimensional image of the colloidal suspension at a 
volume fraction 0.60 acquired by confocal microscopy. The inset 
shows a few typical particle trajectories, from which it is clear that 
we obtain the DOS before large-scale rearrangements happen. 
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FIG. 2: LEFT: Mean square displacement versus time at different 
volume fractions cj> = 0.54,0.57,0.58,0.59,0.60. RIGHT: The 
pressure, as obtained in terms of the volume fraction from the em- 
pirical hard-sphere relation, plotted as a function of the plateau value 
of the mean-square displacement. The full line is a P 3 / 2 fit. 



function of time shows a plateau, indicating that in our exper- 
imental time-window, each particle moves in a 'cage' consti- 
tuted by its neighboring particles. The plateau value decreases 
with increasing volume fraction of the suspension, because the 
particle motion is more restricted (Fig. [2]). Clearly, measure- 
ment times have to be long enough for soft modes (if present) 
to be visible, yet short enough that the system stays within a 
basin. The very existence of a plateau in our measured mean- 
squared displacement shows that there is a range of time scales 
for which local equilibrium can be justified and the soft modes 
can be measured, without many activated events taking place. 

If there is a proliferation of soft modes above a gap, the 
value of the plateau in the mean-squared displacement, as 
Brito and Wyart ifTill showed, is dominated by the modes 
just above the gap. An estimate of the gap value then 
leads lfl4l [l5tl to the prediction that the plateau value of the 
mean-squared displacement scales with pressure as p~ 3 ' 2 . 
This explicit scaling can be tested directly in our experiments 
by plotting the value of the plateau in the mean square dis- 
placement curve vs. the pressure; the latter can be com- 
puted from the volume fraction using an empirical equation 
of state for the colloidal hard spheres, for example llql . 
p = -{k B T/v)(t> 2 ±ln[{{ct> m a X /<t>) {xm - I) 3 ], with v the 



volume per particle. Fitting the data, it is clear that they are in- 
deed compatible with an exponent —3/2 (Fig. [2]). If we admit 
that the system is in local equilibrium during the time-interval 
associated with the plateau, then the values of the squared- 
displacements are independent of the dynamics, which may 
go from overdamped to purely ballistic. (What does depend 
on the nature of the dynamics is the actual time-dependence 
of the displacements, but not their statistical distribution). 

To obtain the Density of States (DOS), we can compute 
the normal modes of the (not necessarily harmonic) vibra- 
tions in the cages from the displacement correlations. De- 
noting u a (t) the components of the particle displacements 
from the average position along the confocal plane u a (t) = 
{(xi(t) — ( xi)) : (y j(t\ — (yi))}> we obtain the displacement 
correlation 1171 ll8L Il9ll matrix (of dimension twice the num- 
ber of observed particles) as, 

D ab = (u a {t)u b (t)), (1) 

where (•} denotes average over the time of the whole period of 
measurement, which is about 175sec. Diagonalizing D ab we 
obtain the normal modes, and the corresponding eigenvalues 
A a . We express the results in terms of the relevant quantities, 

= y/T/K, (2) 

which are the temporal frequencies the system would have ;/ 
it were harmonic. Let us stress that we do not assign any 
meaning of temporal vibrational frequency to uo a at this stage. 

Let us now turn to the results DOS. Because "hard" modes 
are expected to have eigenvalues proportional to the pressure, 
we scale out this trivial effect by plotting the DOS in terms of 
oj/p. The rescaled density of states shows a pronounced peak 
with increased (Fig. [3j> volume fraction in agreement with the 
hard sphere theory [fl4l \vh . 

What do the soft modes look like ? The two dimensional 
eigenvector corresponding to single soft modes shows the 
emergence of rather large vortex-like structures (Fig. |4). On 
the contrary, for higher eigenvalues the eigenvector field looks 
random. To quantify the order in the low-frequency eigen- 
vector field, we calculated the orientational correlation func- 
tion defined as the scalar product (vi.Vj), where Vi and Vj 
are the two-dimensional eigenvector components correspond- 
ing to the i th and j th particle for a single eigenmode. The 
result for eigenmodes in the different parts of the spectrum 
are shown in Fig. [5] Softer modes involve motion correlated 
over many (tens of) interparticle distances, while the hardest 
modes are essentially pairs of particles vibrating with oppo- 
site phases. The correlations characterizing soft modes are 
very weakly dependent on pressure (Fig. |5), suggesting that 
this correlation length stays finite in the large-pressure limit. 

Mode-tomography: In our experimental setup we measure 
the displacements of particles in a small and two-dimensional 
slice of the real system, and we diagonalise a submatrix of the 
full correlation matrix. This is not only due to an experimen- 
tal, but also to a numerical limitation: diagonalizing the full 
correlation matrix is an impossible task. The problem is not 
hopeless, because our slice is typical of all others and will give 
statistically equal results. In particular, the mean squared dis- 
placement per particle and per spatial dimension is the same 
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FIG. 3: Rescaled Density of states at different volume fractions, 
(j> = 0.54,0.57,0.58,0.59,0.60. The frequencies along the hori- 
zontal axis are scaled by the pressure (estimated from the volume 
fraction) times (2cr) 1//3 . The Y-axis is scaled to overlap the high fre- 
quency tails of the distributions. The spectra are consistent with the 
hypothesis that the curves have a low-frequency gap that closes as 
the pressure grows, and approach a universal envelope curve as the 
pressure goes to infinity. 




FIG. 4: Two eigenmodes for volume fraction 0.60. LEFT: low u>: 
collective "swirling" motion is visible. RIGHT: high u motion has 
short range correlations. 
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FIG. 5: An analysis of eigenmodes: the correlation function v\ • vj 
of displacements vl in a single mode (cfr Figure |4](, versus the dis- 
tance 8r = \fl — fj \ between particles. The correlations are com- 
puted for eigenmodes in different parts of the spectrum, as shown in 
D(lo) in the inset by the respective colored lines. Soft modes have a 
long correlation tail, characterizing the "swirling" motion. The hard- 
est modes are made of pairs of particles moving in anticorrelated 
manner. RIGHT: "soft" modes corresponding to systems at differ- 
ent packing fractions. The general shape of the modes seems very 
weakly dependent on the packing fraction. The volume fraction is 
0.60. 

when restricted to a slice and in the whole system. One can 
say more: consider an eigenvector V with eigenvalue A of the 
restricted covariance matrix, involving only the measured par- 
ticles and their motion along the confocal plane. The eigen- 
vectors of the full system {V a } form a complete basis so we 
can expand V in terms of V a as: 

V = J2 c « V« with J2 c * = L (3) 

a a 

One has A = ^ a c„A a , which can also be written as 
J2 a c a ~ ^a) = 0. It is clear from this and the fact that 
eigenvalues are positive that the existence of a large eigen- 
value A for the restricted problem implies the existence of a 
large eigenvalue A Q > A for the complete problem, or, equiv- 
alently, uj a < Q, Each soft mode of the restricted problem 
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implies the existence of - and has a large projection over - 
soft modes of the complete system. In particular, a mode with 
Q = in the restricted problem implies a mode u> — in the 
true system. Although we cannot prove at this point that the 
spectrum associated with a large two-dimensional slice coin- 
cides exactly with that of the full system, the soft modes we 
observe are indeed a reflection of all the soft modes of the full 
system. 

Another argument that justifies the observation on a 
restricted slice is that the vibration we observe in our slice 
is induced by the global soft modes. As mentioned above, 
the DOS D(u>) has a sharp peak at a gap value oj*{p), 
below which there are very few modes (uj*(p) goes to zero 
as p -> oo). The mean squared displacement is domi- 
nated at high pressures by the peak value contribution lfl4ll 
< (8r) 2 >~ D , ' , This means that at large pressure, 
essentially all the motion of each particle is given by a 
combination of the soft modes just above the gap, i.e. those 
having uj ~ uj*{P). 

In conclusion, we have observed an excess of low frequency 



modes in colloidal hard sphere system around the glass tran- 
sition volume fraction. These modes show large-scale cor- 
relations in the velocities of particles, extending over many 
particle diameters. The existence of soft modes had been pre- 
dicted for ideal hard spheres on the basis of theoretical con- 
siderations. Our experimental study shows that they do exist 
colloidal suspensions, probably the most studied system in the 
field of glassy dynamics. These vibrational modes signal the 
onset of macroscopic elasticity and give a microscopic insight 
in the collectivity of the particle dynamics near the glass tran- 
sition. 

We have argued that the normal modes of a restricted set 
of particles gives information on the whole system. Thanks to 
this, a well established technique such as confocal microscopy 
can be used to study global vibrational properties of the sys- 
tem, which contain very detailed information on the geometry 
of the configurations. The above strategy can be applied for 
much more complex colloidal systems. Following particles 
for longer times will further allow us to assess the relevance 
of the low frequency soft modes to long time, activated dy- 
namics near glass transition. 
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